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The zero recoil limit of the B — ► Dlv form factors is calculated on the lattice, which provides a model- 
independent determination of |V c f,|. Considering a ratio of form factors, in which the bulk of statistical and 
systematic errors cancel, we obtain a precise result both for h+(l) and for h-(l). 



1. Introduction 

For the determination of \V c b\ through the ex- 
clusive decay B — > D^*Hv, the theoretical calcula- 
tion of the form factor T B ^ D (»)(w), especially in 
the zero recoil limit, is necessary. Previously this 
has been done using the zero-recoil sum rule [Q 
or the heavy quark expansion of the form factor 
Hj^|. Both of these calculations, however, need to 
introduce an assumption or a theoretical model 
to deal with hadronic effects away from the in- 
finite quark mass limit. In this talk, we present 
a lattice calculation of the form factor Fb^dO-), 
which can be used for a model independent de- 
termination of | V c b | , including deviations from the 
heavy quark limit. 



2. B — > Dlv Form Factors 

The differential decay rate of B 
portional to the square of 



Dlv is pro- 



T b ^d{w) = h + (w) 



m B - mp 
tub + mo 



h-(w). 



(1) 



h+(w) and h-(w) are form factors defined 
through 

(£>(p')|Vl£(p)> = Vm B m D 

x [h+(w)(v + v'Y + h-(w)(v - v'Y], (2) 

where — p^/ms, v'^ = p'/rriD and w = v ■ v' . 
In the heavy quark mass limit, h-(w) vanishes, 
and h+(w) agrees with the universal form factor 
£(u>) (the Isgur-Wise function), which is normal- 
ized in the zero recoil limit £(1) = 1 0]. 



The I/rag expansion may be used to describe 
the deviation from the heavy quark limit. At zero 
recoil, the expansion becomes 
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h-(l) = 0-c 



c w (— 

















0(l/m 3 Q ),(3) 



+ 0(l/m 2 Q ). (4) 



The absence of the O^I/ttiq) term in the expan- 
sion of h+(l), which is a part of the Luke's theo- 
rem JHJ , is a consequence of the symmetry under 
the exchange of m c and m& (see eq.(^)), and this 
particular form of the expansion is also restricted 
by the symmetry (anti-symmetry for /i_(l)). Our 
task is to determine the coefficients Cj? and c_ , 
for which there has been no model independent 
calculation. 

3. Mi) 

The vector current V M = cj^b appearing in 
eq.(||) must be related to the lattice counter- 
part vl att using the perturbative relation = 
ZyVjf tt '. This is true even for the equal mass 
(m c = nib] case, because the lattice (local) vec- 
tor current is not conserved. Without a two-loop 
calculation, this perturbative matching could in 
principle be a source of large systematic uncer- 
tainty of O(cVj) ~ 5%, which is too large to ob- 
tain the precision we seek for the form factor 
(< 5%). The statistical error in the lattice cal- 
culation would also be a problem, if we employed 
the usual method to extract the matrix element 
from three-point correlator. 
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Figure 1. R (t) as a function of t. 

To reduce these errors, we define a ratio at zero 
recoil 

-i cont 



R 



B^D 



(D\V \B)(B\V \D) 



l(D\Vo\D)(B\V \B) 
D {l)h^ B {l) 



D (l)h B - B (l) - 
used the property 



\h» 



(5) 



where we 

h B ^ B (l) = 1 derived from current conservation. 
R B ^ D may be related to the lattice counterpart 



Z~y cb Z-ybc 

Zv~Z- 



ybh 



(D\Vo\B)(B\V Q \D) 



-\ latt 



D\V \D)(B\Vo\B)\ ' (6) 

The ratio of matching factor can be safely evalu- 
ated with perturbation theory, since a large can- 
cellation of perturbative coefficients takes place in 
the ratio. The one-loop calculation is discussed in 
a separate talk Q . 

To calculate the hadronic amplitude, we define 

c DV aB{t)c BV D {t) 



R 



(*) = 



R 



B^D 



(7) 



C DV ° D (t)C BV ° B {t) 

where C DV ° B (t) is a three-point correlator, whose 
initial and final states are fixed at t = and t = 
Nt/2 respectively and the vector current is moved 
to find a plateau. Our calculation is done on a 
12 3 x24 lattice at (3 = 5.7. The Fermilab action 
is used for heavy quark with c sw — 1/uq. 

Figure |l| shows a nice plateau for the ratio 
R B ^ D (t) . Even with only 200 configurations the 
statistical error is remarkably small (< 1%), be- 
cause of the cancellation of statistical fluctuations 
in the ratio. 

Fitting the plot with a constant we obtain 
|/i+(l)| 2 for each combination of initial and fi- 
nal heavy quark masses. To fix the parameter in 
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Figure 2. (1 - \h+(l)\)/(l/m c - 1/rrib) 2 versus 
l/am c + l/arnb- 

the 1 / tuq expansion, we choose six values of the 
heavy quark mass covering the physical rrib and 
m c , and fit the results with the form 



h+(l) 
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(8) 



where the 0(1/toq) term is required to explain 
the data. Figure shows the quantity (1 — 



ft+(l))/(l/m c - 1/nib)' 2 = cf - c^(l/m c + 
1/irib), from which we extract the coefficients c + 

= 0.029(11) and 
= 0.011(4). In physical units we obtain 
cy = (0.20(4)GeV) 2 and c { + ] = (0.26(3)GeV) 3 . 
Shifman et al. [Q derived a bound > — 
P? G )j2 = (0.26±^GeV) 2 using the zero-recoil 
sum rule. ^ Our result is consistent with this 
bound within errors. 
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4. mi) 

h-(l) cannot be obtained from the matrix el- 
ement at zero recoil. We introduce a finite (but 
small) D meson momentum p' and define a ratio 

B ^ D = (Dip'MlBjO)) 
v </ y ° (D(p')\V \B(0)) 
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h-(w) 
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h-(w) 
h+(w) 



(9) 



1 For the hadronic parameters we used ^ = 0.5(l)GeV 2 



and 



0.36GeV 2 
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Figure 3. 



R 



t 



(i) for the final 



state momentum (1,0,0) (circles) and (2,0,0) 
(squares). B —> D (filled symbols) and D —> B 
(open symbols) are shown. 



where we use the definition of the form factors 
and expand in small v' 2 . It is more convenient to 



define a double ratio 



R 



(B^D)/(D^D) 
Vi/V 



U V t /V l n Vi/V 



h-(w) 
h + (w) 



1 h-(w) ,2 
1h+(w) V 



(10) 



The property of elastic scattering h D \~^ D (w) = 
is used for the denominator. 

Figure shows the corresponding ratio 



R 



v .iy (t) on the lattice for two differ- 

ent values of the D meson momentum. As in the 
calculation of /i+(l), the plateau is very clear and 
we can extract (l — h—(w)/h+(w)) from this plot. 
The small correction proportional to v' may 
be eliminated by extrapolating to the v' 2 — * 
limit. We also observe an antisymmetric prop- 
erty h D ^ B (w) = -h B ^ D (w). 

The heavy quark mass dependence can be ob- 
tained with a similar strategy. We fit our data 
with 
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and obtain c 



(i) 



(2) 



0.23(3) and c [ _ 

(i) 



(11) 



0.06(1), 



which correspond to cL — 0.26 (4) GeV and c 
(0.29(3)GeV) 2 in physical units. 



J2) 



5. F b ^d{1) 

The results with physical mass parameter are 

Ml) = 1-016 ±0.003 ±0.002 ± 0.006, (12) 
Ml) = -0.112 ± 0.014 ± 0.011 ±0.025, (13) 

where h + (l) includes the one-loop correction 
±0.025(6) ^). Errors arise from statistics, mass 
parameter determination, and our estimate for 
higher order perturbative correction in the order 
given. Using (0) we obtain 

Fb^d{±) = 1-069 ± 0.008 ± 0.002 ± 0.025. (14) 

6. Conclusions 

Using a ratio in which a large cancellation of 
statistical and systematic errors takes place, we 
have calculated the B — > Dlv form factors very 
precisely, which may lead to a better determina- 
tion of \Vcb\- In the error bar given above, how- 
ever, we have not yet included the discretization 
error and the effect of quenching. It is our hope 
that the bulk of these errors also cancels in the 
ratio. We leave these important issues for future 
study. 

Our method can also be applied to the B — * 
D*lv form factor, which currently yields an ex- 
perimentally more precise determination of \V c b\- 
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